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Abstract. Le Brun et al. (1997) presented the first identi- 
fications of the galaxies giving rise to 7 intermediate red- 
shift damped Lya (DLA) absorption systems. Here, we 
study the gravitational lensing properties of these fore- 
ground galaxies based on their observed optical appear- 
ance and on the absence of any secondary lensed quasar 
image. We consider the possibility that any secondary im- 
age be hidden due to extinction by dust, but find it un- 
likely. We derive upper limits on the amplification factor 
affecting the luminosity of the background quasars; in each 
case, this factor is found to be less than 0.3 mag. We also 
obtain upper limits on the total mass of the damped Ly-a 
galaxies, within radii equal to the quasar impact parame- 
ters. Mass-to-light ratios are found to be consistent with 
existing estimates based on X-ray emission or on motion 
of dwarf satellites. Although we show that lensing is not 
important in this sample, we note that existing DLA sur- 
veys used to determine the cosmological density of gas at 
z < 1 are based on samples of quasars brighter than the 
ones considered here and for which the amplification bias 
is likely to be stronger. 

Key words: Gravitational lensing - Galaxies: halos - 
Galaxies: ISM - Quasars: absorption lines 



1. Introduction 

Redshifted damped hydrogen Lya (DLA) lines are some- 
times detected in the spectra of background quasars. 
They present neutral hydrogen column densities A^(H i) > 
2 x 10 20 cm~ 2 similar to the value found at the optical 
radius of present-day galactic disks (Wolfe et al. 1986), 
which suggests that they arise in progenitors of present- 
day spiral galaxies. 
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Because the impact parameters necessary to produce 
such absorption lines arc likely to be small, the galaxy 
gravitational potential may lead to significant gravita- 
tional lensing effects, such as amplification of the apparent 
luminosity of the background quasar, or even formation of 
multiple images. 

In turn, these effects can induce severe biases in the 
samples of quasars and absorbers. First, impact parame- 
ters are increased because of lensing (the by-pass effect; 
see Smette, Claeskens, Surdej 1997, hereafter, SCS) thus 
decreasing the cross-section for DLA absorption. Second, 
background quasars which present DLA absorptions in 
their spectra also have their apparent luminosity increased 
because of lensing: hence, intrinsically faint quasars can 
enter apparent magnitude limited samples. As a conse- 
quence, these samples may show an over-representation of 
quasars that present DLA absorption systems: this am- 
plification bias apparently increases the cross-section for 
damped Ly-a absorption. These aspects can be evaluated 
in a statistical sense for a well-defined magnitude limited 
sample (cf. SCS; Bartelmann &; Loeb 1996). For bright 
quasar samples, the amplification bias was found to domi- 
nate the by-pass effect, which leads to predictions that the 
number density dJV/dz of z ~ 0.5 DLA systems and the 
cosmological density of neutral hydrogen f2jj j (since DLA 
absorbers contain most of the H I of the Universe) may be 
significantly overestimated. Note that extinction by dust 
introduces other biases competing with the ones just de- 
scribed (Fall & Pei 1993, Perna et al. 1997, Bartelmann 
& Loeb 1998). 

It must be noted however, that existing theoretical es- 
timates are based on a series of simplifying assumptions, 
one of which being that the DLA absorbers are very sim- 
ilar to present-day disk galaxies. As a matter of fact, Le 
Brun et al. (1997) presented the first identifications of 
damped Lya absorbing galaxies at intermediate redshifts: 
in a sample of 7 DLA absorbers, 3 turn out to be spi- 
ral galaxies of various sizes and luminosities, 2 are com- 
pact objects and 2 are amorphous, low-surface brightness 
galaxies (an additional compact object appears to be re- 



sponsible for a higher-redshift DLA system). This study 
showed that the population of galaxies giving rise to DLA 
systems is more heterogeneous than previously thought. 

Independently, SCS found that lensing effects may be 
present in the IUE (Lanzetta et al. 1995) and Rao et al. 
(1995) Mgn surveys. In addition, in a follow-up of this 
latter survey, Rao & Turnshck (2000) find a surprisingly 
high value of Sljjj at z < 1, as expected if lensing effects 
are important (SCS). 

Here we consider the lensing properties of the absorb- 
ing galaxies detected by Le Brun et al. (1997). These ob- 
servations were made using the WFPC2 camera onboard 
the Hubble Space Telescope, so that the absorbing galaxy 
candidates are clearly separated from the quasar images, 
and their characteristics (impact parameter, luminosity, 
morphology) are fairly well estimated. 

The sample of quasars described here is not complete 
in any sense except that it contains all the DLA systems 
known at the time of writing the HST proposal (cycle 6). 
For illustration purpose only, we have evaluated the effects 
of lensing following SCS (exponential model plus Gaussian 
hole, or E+GH). These theoretical estimates are based on 
the observed magnitude 777,450 and emission redshift z e of 
the quasar and on the DLA absorber redshift Z&. Would 
the DLA sample be representative of a well-defined, com- 
plete magnitude limited sample, we find that the number 
density of DLA systems dAf/dz would have been overesti- 
mated by 18%. In other words, statistically speaking, one 
of the quasars in this sample would not have been observed 
if gravitational lensing was not taking place. In addition, 
the value of Sljj l would have been over-estimated by 28%. 
We also find that the probability that one quasar in the 
sample presents multiple-imaging would be close to 50%. 
The values H = 100 km s" 1 Mpc -1 and (0, A) = (0,0) 
are used throughout this paper. 

2. Lensing properties of the DLA absorbing 
galaxy/quasar pairs 

The observed characteristics of the seven quasar/galaxy 
pairs are presented in Table [I]. We present three different 
methods to evaluate the lensing properties of these seven 
configurations: (i) a model-independent formalism; (ii) a 
formalism which assumes a Singular Isothermal Sphere 
(SIS) lens model, only based on the geometry of the sys- 
tems (i.e. redshifts of the quasar and the galaxy in each 
association and their angular separation); (Hi) a formal- 
ism which assumes a SIS model whose velocity dispersion 
is determined from the observed luminosity of the DLA 
galaxies and the Tully-Fisher relation. 

The values of the angular impact parameter B\ and 
absolute luminosity Lb of the absorbing galaxy image were 
obtained by Le Brun et al. (1997) and Boisse et al. (1998). 

No secondary lensed QSO image is detected in any of 
the systems presented here at more than 0.3" from the 
quasar, down to a limiting magnitude of m^ = 25.5, i.e. 



6.8 to 9 magnitudes fainter than the main observed QSO 
image. The absence of any secondary image either means 
that the lensing configuration is indeed not capable to 
produce multiple imaging, or that the apparent luminosity 
of the secondary image is too faint to be detected, possibly 
due to extinction by dust. 

2.1. (Nearly) model-independent constraints 

Subramanian & Cowling (1986) showed that, for a spheri- 
cal mass distribution, with surface density decreasing from 
the center to the outer parts, as expected for individual 
galaxy halos, a sufficient and necessary condition to have 
multiple images is that the central value S is larger than 
a critical value S cr ; t defined by : 



D n 



4ttG L>ol As 



(1) 



and thus independently of the model chosen for the halo. 
In this relation, c is the speed of light, G the gravitational 
constant, D os , D a \ and D\ s the angular-diameter distances 
between the observer and the source (QSO), the observer 
and the lens, the lens and the source, respectively. Further- 
more, they show that the absence of a secondary lensed 
QSO image ensures that So < S cr it still for a spherical 
mass distribution. They conjecture that this result is also 
valid for centrally peaked elliptical lenses. We therefore 
conclude that the absence of a secondary QSO image at 
a separation larger than 0.3" in our sample ensures that 
S < S C rit over such an angular scale. 

On the other hand, for a lens with circular symmetry, 
the mean surface density within the Einstein radius Be 
is equal to E cr it (cf. Schneider, Falco & Ehlers 1992); in 
addition, the location of the main image is always such 
that its angular separation from the lensing galaxy 6\ > 
#e- In particular, this relation is true even in the case of 
multiple images with any secondary image being hidden 
due to extinction by dust. Consequently, we can derive an 
upper limit on the projected mass M,^ ax (< n) enclosed in 
a disk of radius r\ = D \ 9\ centered on the galaxy, as 
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and thus also on the average M/L ratio within rj. 

No constraint can be set on the amplification A of the 
QSO image. 

2.2. SIS lens model with only geometrical constraints 

We here assume that the distribution of matter within 
each galaxy can be described as a singular isothermal 
sphere (SIS), whose volume mass density p is given as 
a function of the distance r to the galaxy center by: 
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Table 1. Characteristics of the seven galaxy/quasar pairs: quasar emission redshift z e , apparent magnitude in the 
F450 filter (similar to the B band), DLA line redshift z^, neutral hydrogen column density iV(Hi), angular impact 
parameter &j, linear projected impact parameter n, fc-corrected absolute magnitude Mb and its estimated la error 
of the absorbing galaxy (see Boisse et al. 1998 for details). The last column provides the galaxy morphological type, 
derived from the WFPC2 images as precisely as possible. 



Quasar name 


Ze 


rri450 


Zd 


JV(Hi) 


0i 

(") 


n 

(h- 1 kpc) 


Mb 


Type 


EX 0302-223 


1.400 


16.5 


1.0095 


20.39 


1.14 


6.2 + 0.1 


-18.1 + 0.1 


Compact Dwarf 


PKS 0454+039 


1.345 


17.6 


0.8596 


20.69 


0.80 


4.1 + 0.1 


-18.5 + 0.1 


Compact Dwarf 


3C 196 


0.871 


18.7 


0.4370 


20.8 


1.10 


4.1 + 0.1 


-20.0 + 0.05 


Face-on Giant Spiral 


Q 1209+107 


2.191 


18.6 


0.6295 


20.2 


1.62 


7.3 + 0.1 


-19.7 + 0.05 


Edge-on Spiral 


PKS 1229-021 


1.038 


17.6 


0.3950 


20.75 


1.40 


4.9 + 0.1 


-16.9 + 0.2 


Diffuse Dwarf 


3C 286 


0.849 


18.0 


0.6920 


21.19 


0.90 


4.3 + 0.1 


-18.7 + 0.2 


Irregular 


MC 1331+170 


2.084 


16.6 


0.7446 


— 


3.86 


18.8 + 0.0 


-20.9 + 0.2 


Edge-on Bright Spiral 



where a v is the 1-dimensional velocity dispersion of the 
SIS. As a consequence, the total projected mass enclosed 
in a disk of radius r is 






(4) 



We then compute the angular Einstein radius for each 
quasar/galaxy configuration, by 
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In the case of an SIS, the observed images are located at: 



9i = 9 s ± 



(G) 



where 9s is the true (unobserved) position of the source. 

Multiple (double) imaging only occurs for 0\ < 26s, 
in which case the separation between the two images is 
A6> = 20 E , and their amplification is Ai = {9e/0s) ± 1. If 
6\ > 2#e only one image is formed and its amplification is 
Aj = (0 E /0 S ) + 1. 

Inversely, even in the case of an isothermal sphere with 
a core-radius, Narayan & Schneider (1990) showed that 
the condition 9\ > 2#e is sufficient to avoid the formation 
of multiple images. Using only the geometry of the system, 
we can set the following constraint on the mass of the 
DLA galaxy. We successively consider the cases of single 
and double image systems : 

2.2.1. Single image system 

If only one image is present, 9\ > 2 #e- Using Eq. and 
inverting Eq. imply that: 



Mgcom / . \ ^ 
flS.l (< r l) < T^ 
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The only constraint on the amplification affecting the ob- 
served image is A < 2. 



2.2.2. Double image system 

Suppose now that two images are actually produced by 
the lens. In order to explain the observations, we have to 
assume that only the brightest image is detected, while the 
faintest one is affected by dust extinction by an amount 
Eb so that its apparent magnitude is fainter than the lim- 
iting magnitude of the corresponding WFPC2 observa- 
tions. 

The angular separation between the observed image 
and the center of the lensing galaxy B\ is such that #e < 
6\ < 2#e- Consequently, the total projected mass located 
within a disk of radius r\ is constrained by : 



C(< rO 
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Let's define Am b s as the difference between the 
limiting magnitude of the corresponding WFPC2 frame 
and the magnitude of the main observed QSO image : 
Am bs = "i45o — "i45o(QSO). We then have the following 
relations between the amplifications of the two images I 
and B of the quasar : 



A 
~Ai 



£2. = 15. < 10 -0.4[Am obs -(E B -Bi)] 



(9) 



where E\ and £b are the extinction affecting the main 
and secondary images, respectively. 

Since the secondary image can be very close to the 
center of the deflecting galaxy, it might be extremely ex- 
tinguished, and thus undectable on our data. 

2.2.3. Summary 

Only using geometrical quantities with a SIS lens model 
does not greatly improve the limits based on a model- 
independent formalism: the value M^ ax (< n) is within a 
factor of two the upper limit of the mass of a SIS lens 
only constrained by the angular separation between the 
lensing galaxy and the quasar, independently of the fact 
that multiple images exist or not. 



2.3. SIS lens model constrained by the Tully-Fisher 
relation 

In this section, we consider one additional piece of infor- 
mation brought by the HST observations: the luminosity 
of the galaxies considered to be responsible for the DLA 
absorption. 

If we assume that the 1-dimensional velocity disper- 
sion of the SIS cr v is equal to the velocity dispersion of 
matter in the galaxy a^ F , which is related to the maxi- 
mal value i> m ax of the rotational velocity of a galaxy by 
crJ F = V ma x/v2j we can indeed use the Tully-Fisher re- 
lation (Tully & Fisher 1977) to derive its value from the 
galaxy luminosity: 



TF 
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(10) 



where cr*' TF = 144±* 3 km s~\ M±(B^) = 
— 19.9^0 4 + 51og/i (h is the Hubble constant in units of 
100 km s _1 Mpc -1 ), and «tf = 2.6±0.2; these values are 
in accordance with the work by Fukugita & Turner (1991) 
and SCS. This relation is derived from local galaxies, but 
most recent studies (see Bershady 1996 for a review) in- 
dicate that there is no evolution of the parameters in this 
relation at intermediate redshifts. Then, even if the value 
of M* is known to evolve with redshift (see e.g Lilly et 
al. 1995), this only indicates an evolution of the density 
of galaxies of a given luminosity, but not of the relation 
between the luminosity and the velocity dispersion of a 
given galaxy. 

We can then set an upper limit to the projected mass 
enclosed in a disk of radius n centered on the galaxy, as 
Eq. becomes 



M s \U<n) = 



* K TF ) 2 
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We can also estimate the amplification factor Aj F of 
the quasar apparent luminosity due to gravitational lens- 
ing by the galaxy responsible for the DLA absorption, as 
0g F is derived by inserting o~^ F in Eq. |5j: 
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The values of Afgj F (< rj) and Aj F derived for each 
of the QSO-galaxy association are presented in Table 

2. As expected, the inferred mass-to-light ratios Mgj F (< 
r^/L^, also listed in Table 2, are close to the mass-to- 
light ratios for spiral and elliptical galaxies at the present 
epoch estimated by various, classical methods (see Bah- 
call et al 1997 for a review). This result, discussed further 
in Section 3.3, ensures us that the values obtained for the 
amplification factor Aj F are probably good estimates. 

3. Discussion 



3.1. Absence of a secondary image 

As can be seen in Table 2, application of the Tully-Fisher 
relation to the absolute luminosity of the DLA galaxies 
implies that the impact parameter of the line-of-sight to 
the QSO never falls within twice the value of the galaxy 
Einstein radius. Hence, we do not a priori expect any sec- 
ondary lensed QSO images, which is confirmed by the 
HST/WFPC2 observations. 

However, the total luminosity of the galaxies might be 
underestimated, because of extinction due to the presence 
of diffuse dust in the galaxies themselves (self-extinction), 
especially in the rest-frame B band in which our galaxies 
are observed. Since the Tully-Fisher relation given above 
in Eq. 10 was determined from a sample of local galax- 



ies whose B magnitudes were corrected for self-extinction, 
even for inclined systems, we are led to use a global ex- 
tinction. 

In particular, a significant self-extinction correction 
might be applied to a DLA galaxy absolute B magnitude 
if it presents a large inclination. Thus, from the observed 
values, we can calculate the self-extinction which is neces- 
sary to give a real Einstein radius large enough to lead to 
Q\ < 26*e and the formation of a double image. For four ab- 
sorbing galaxies, this value of the self-extinction is larger 
than 3 magnitudes, which we consider to be implausible, 
given that such values have not been observed in samples 
dedicated to the study of self-extinction in galaxies (Xu 
et al. 1997). 

For the three remaining galaxies (toward 3C 196, Q 
1209+107 and MC 1331+170), the self-extinction neces- 
sary to lead to a value of 8e compatible with multiple 
imaging is smaller, ranging from 1 magnitude for 3C 196 
to 1.6 and 1.8 magnitude for the last two. The last two 
values are comparable to the highest ones detected in the 
sample of Xu et al. (1997). On the other hand, we con- 
sider unlikely that the galaxy responsible for the DLA 
in the spectrum of 3C 196 is affected by 1 magnitude of 
self-extinction: the galaxy is seen face-on and its redshift 
is z d = 0.4370, so that the F702W band is actually cen- 
tered at 4900A in the galaxy rest-frame. Furthermore, the 
galaxy extent as determined in the F450W image, which 
roughly corresponds to near UV in the galaxy rest frame, 
is nearly equal to the one seen in the F702W filter. The 
self-extinction correction is thus probably small in this 
object. 

In summary, the self-extinction required to produce 
multiple imaging in these 3 systems are unlikely to be sig- 
nificant. However, they are not unrealistic. Let's assume, 
therefore, that multiple imaging is taking place. Conse- 
quently, the absence of detected secondary images allows 
us to provide some constraints on the extinction in these 
galaxies along the particular lines-of-sight to the QSO im- 
ages. 

Indeed, if multiple imaging is taking place for 3C 196, 
Q1209+107 and MC1331+170 and if the SIS model is an 
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Table 2. Lensing properties of the seven quasar/galaxy pairs. Symbols are defined in the text. The TF upper-script 
indicates that the corresponding quantity has been evaluated with the Tully-Fisher relation. Calculations have been 
done using Hq = 100 km/s/Mpc, il = 0, A = 0, and no correction for extinction. 
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adequate representation of the matter distribution within 
the lens, then the observations could reveal two images 
with a magnitude difference Atobi smaller than Am b s ,j — 
AE, where i = 1, 2, 3 represents each of the three quasars. 
From Eq. 9, such a situation occurs if 



#s < fe &e, 
where 



fo = 



1 -j. lQ-°- 4 Am BI ' 



(13) 



(14) 



The probability Pi(AE; Am b s ,i) of not detecting a 
secondary image in the system i is thus given by 



Pi = l- 



y E 



1 - (fe? 



(15) 



This last relation assumes that the sources are uni- 
formly distributed behind the lenses: we neglect the am- 
plification bias, which tends to select systems with small 
#s and, consequently, to strengthen the following conclu- 
sion. 

The probability P{AE; Am obs ,i, Am obs ,2, Arn obs ,3) of 
not detecting a secondary image in any of the 3 systenrs 
is therefore 



4=1 

which is represented as a solid line in Fig. |l|. 



(16) 



We can see that the non-detection hypothesis, i.e. the 
observations, is ruled out with a confidence level larger 
than 3 sigma if the differential extinction is smaller than 
3.9 magnitudes on each sightline. But an extinction larger 
than 3.9 mag is only expected in very dense clouds, whose 
covering factor is very small. 

Furthermore, Table 2 shows that the limits on the 
mass-to-light ratios inferred from lensing and geometrical 
constraints alone are usually an order of magnitude larger 
than the ones derived from an application of the Tully- 
Fisher relation. If a secondary image was hidden due to 
extinction, it would imply that the mean M/L ratio of 
galaxies is much higher than previously thought. 

For all these reasons, we maintain our conclusion that 
the observed configurations of DLA absorbers are not 
likely to give rise to multiple images. 



3.2. Biases due to gravitational lensing 

A quasar located at an impact parameter slightly larger 
than the Einstein radius can still have its apparent lumi- 
nosity increased by a factor A, due to gravitational lensing 
amplification. Column 10 of Table presents this factor 
estimated for each of our systems: as one can see, it is al- 
ways smaller than 0.3 mag., with an average value of 0.12 
mag. The values for a (O, A) = (0.3, 0.7) Universe are only 
slightly larger (by 3 to 6%). 



-2 



-4 




3.3. Comparison of the M/L ratio constraints with other 
methods 



AK (mag) 

Fig. 1. Probability P of not detecting any secondary im- 
age for the three selected quasars (see text) as a function of 
the differential extinction AE toward the three sight- lines. 
The horizontal dashed line represents the 2>a (or 0.25%) 
level, which we choose as a threshold 



As we now have an estimated value for the amplifi- 
cation, instead of estimating statistical lensing effects, we 
can actually compute lensing effects for each QSO individ- 
ually. We used the approach of Narayan (1989) to evalu- 
ate the excess of quasars close to foreground galaxies. This 
method takes into account both the amplification bias and 
the by-pass effect. We point out, however, that it aims at 
estimating the excess number of quasars in the vicinity 
of galaxies, not the excess of quasars in the vicinity of 
galaxies giving rise to DLA systems in quasar absorption 
spectra. However, we prefer this method due to uncertain- 
ties in the determination of the inclination of the galaxies 
and other observational variables, and because it is good 
enough for our purpose. 

The mean excess of quasars close to the DLA galaxies 
is found to be equal to a factor ~ 1.14 for a (fi, A) = (0, 0) 
Universe. Therefore, if these quasars were drawn out of a 
magnitude limited sample, 14% of the quasars that present 
a DLA line in their spectrum would have been observed 
because of gravitational lensing (this calculation should 
be considered as a mere exercise, as the sample is not 
complete). This value is close to the 18% obtained by the 
method described in SCS (cf. Sect. 1). 

On the other hand, if these quasars were drawn out of 
a volume limited sample, the magnification bias would be 
irrelevant: only the by-pass effect would be acting. In this 
case, the fact that the observed impact parameters 9\ are 
larger than 9a, leads to an underestimate of f2jj j by about 
6%, based on the E+GH model of SCS. 



We have plotted in Fig. the upper limits to the model- 
independant mass-to- light ratios M,„ ax (< r^/Lg* we de- 
rive for the seven galaxies. The dotted and dashed lines 
show the power laws representing the overall mass-to- light 
ratio in spiral and elliptical galaxies, respectively, as esti- 
mated from Bahcall et al. (1995): this review presents re- 
cent results obtained on the mass-to-light ratio using the 
classical methods (H I emission, X-ray emission, motion of 
dwarf satellites ...). 

As can be seen, the strongest limits set by gravita- 
tional lensing, which are obtained for the three spirals of 
the sample, are an order of magnitude above the values 
obtained from other methods. Thus, although they are 
compatible with other estimations, these observations do 
not allow to efficiently constrain the value of the hidden 
mass in the galaxies responsible for the DLAs. 

3.4. Future work 

We note that the published surveys to determine the cos- 
mological density of neutral hydrogen Sljj l at < z < 1 
have been carried on using samples of quasars that are 
generally brighter than all the quasars on which the 
present study is based: only 15% of the quasars in the 
IUE survey (Lanzetta et al. 1995) have a B magnitude 
fainter than the quasars presented here; half of the Rao et 
al. (1995) quasars (followed up by Rao & Turnshek (2000)) 
are brighter than the brightest quasar in our sample, and 
none are fainter than the faintest quasar in this same sam- 
ple (some quasars are actually used both in this paper and 
in Rao et al.'s). These surveys may thus be more strongly 
affected by the magnification bias than the sample pre- 
sented here; consequently, the results of this paper should 
not be interpreted as meaning that lensing effects have 
negligible effects on surveys of DLAs at < z < 1 . 

In order to have a better constraint on the lensing ef- 
fects in DLA surveys, we have carried a HST-NICMOS 
survey of 13 bright quasars whose spectra present a DLA 
system at low redshift, including the sample presented in 
this paper. 
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